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Abstract 

Canonical ensembles consisting of M-unit Hubbard dinners have been studies within the nonex- 
tensive statistics (NES). The temperature dependences of the energy, entropy, specific heat and 
susceptibility have been calculated for the number of dimers, M = 1,2,3 and oo. We have assumed 
the relation between the entropic index q and the cluster size given hy q = 1 + 2/N (N = 2 M 
for M dimers), which was previously derived by several methods. For relating the physical tem- 
perature T to the Lagrange multiplier /3, two methods have been adopted: T = l/ksP in the 
method A [Tsallis et al. Physica A 261, 534 (1998)], and T = Cq/ksP in the method B [Abe et 
al. Phys. Lett. A 281, 126 (2001)], where ks denotes the Boltzman constant, Cq = YliPi^ 
Pi the probability distribution of the ith. state. The susceptibility and specific heat of spin dimers 
{Heisenberg dimers) described by the Heisenberg model have been discussed also by using the NES 
with the methods A and B. A comparison between the two methods suggests that the method B 
may be more reasonable than the method A for nonextensive systems. 
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I. INTRODUCTION 



In the last several years, much study has been made with the use of nonextensive statistics 
(NES) which was initiated by Tsallis [l-sl. When the physical quantity Q of a system 
consisting of particles is expressed by Q oc A^"^, it is called intensive for 7 = 0, extensive 
for 7 = 1, and nonextensive for 7 7^ 1 and 7 7^ 0. For example, in a spatially homogenous d- 
dimensional classical gas with the attractive interaction decaying as r~", we get 7 = 2 — a/d 
for < a/c? < 1 (nonextensive) and 7 = 1 for a/d> 1 (extensive) [3]. The nonextensivity is 
generally realized when the range of interactions is long enough compared to the linear size 
of the system. Then small-scale systems may be nonextensive even when the interaction is 
not long-ranged one. 

Tsallis has proposed the NES entropy given by [ll 

S, = kB , (1) 



l-q 

where fc^ is the Boltzman constant, q the entropic index, and pi the probability density of 
the ^th state. Note that the entropy of BGS is obtained from Eq. (1) in the limit of g = 1. 
The quantity | g — 1 | expresses the measure of the nonextensivity. The NES has been 
successfully applied to a wide range of nonextensive systems including physics, chemistry, 
mathematics, astronomy, geophysics, biology, medicine, economics, engineering, linguistics, 
and others 

In our previous papers , we have applied the NES to Hubbard dimers described by 
the two-site Hubbard model. The Hubbard model is one of the most important models 
in solid-state physics (for a recent review, see Ref. [71). The Hubbard model consists of 
the tight-binding term expressing electron hoppings and the short-range interaction term 
between two electrons with opposite spins. The Hubbard model provides us with good qual- 
itative description for many interesting phenomena such as magnetism, electron correlation, 
and superconductivity. In particular, the Hubbard model has been widely employed for a 
study on transition-metal magnetism. Thermodynamical properties of grand-canonical en- 



sembles of a single Hubbard dimer have been calculated within the NES jsj. It has been 
shown that specific heat and susceptibility calculated by the NES may be significantly dif- 
ferent from those calculated by the Boltzman-Gibbs statistics (BGS) when the entropic 
index q departs from unity, the NES with q = 1 reducing to the BGS. It is interesting to 
compare the calculated results with experimental data. However, experimental data for a 



single dimer as adopted in Ref. jsl, is not available in actual experiments. Usually exper- 
iments on nanosystems are performed on samples which include many clusters consisting 
of, for example, multiples dimers (for reviews, see Refs. 8-10). Iron S = 5/2 dimers (Fe2) 
in [Fe(OMe)(rf6m)2]2 



ll| have the nonmagnetic, singlet ground state and their thermody- 



namical property has been analyzed with the use of the Heisenberg model 12l| - ll4l . Similar 
analysis has been made for transition-metal dimers of V2 [iS^, Cr2 [iG^, Co2 [o], Ni2 [l^ 
and Cu2 [l9^. Some charge-transfer salts like tetracyanoquinodimethan (TCNQ) with dimer- 
ized structures, have been analyzed by using the two-site Hubbard model within the BGS 



20| . Their susceptibility and specific heat were studied by taking into account the inter- 



dimer hopping, whose effect is negligibly small 2C|. Such procedure may be justified within 
the BGS where the specific heat and susceptibility are treated as the extensive quantities: 
macroscopic measurements are expected to reflect the property of a constituting dimer. This 
is, however, not the case in the NES. 

The purpose of the present paper is two folds. 
(1) It is interesting and indispensable to investigate how thermodynamical properties may 
change when the size of a given cluster is varied within the NES. It has been shown by 
several methods that the entropic g of a nanosystem consisting of independent particles 



is given by 



2l|-|23| 



q = l + 



N' 



(2) 



Bearing in mind a magnetic cluster containing M transition-metal dimers, we have employed 
the Hubbard model to perform NES calculations for various M, assuming the relation given 
by 



q = l + 



M' 



(3) 



which is derived from Eq. (2) with = 2M for M dimers. We have adopted 1.0 < g < 2.0 
for 1 < M < oo in this paper, preliminary results of the M dependence of thermodynamical 
quantities of clusters having been reported in Ref. gJ. 

(2) It is not clear in the current NES how to relate the physical temperature T to the 



Lagrange multiplier /3 



51]. The following two methods have been so far proposed: 
1 



T 



kBf3' 



(method A) 
(method B) 



(4) 
(5) 
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where Cg = ^jP^ The method A proposed in Ref. 2| is the same as the BGS. The 
method B is introduced so as to satisfy the zeroth law of thermodynamical principles and 



the generalized Legendre transformations 2J]. It has been demonstrated that the negative 



specific heat of a classical gas model which is realized in the method A [25], is remedied in 



the method B 



24| . The specific heat and susceptibility calculated by the two methods A 



and B are qualitatively similar, but quantitatively different: the nonextensivity calculated 
by the method A is generally more significant than that calculated by the method B js, 6|. 
In particular, the Curie constant of the Hubbard model in the limit of vanishing couplings 
calculated by the method A becomes spuriously large {sgI while that calculated by the 
method B is reasonable jsl- This is consistent with the result for localized free spins 5]. In 
order to get more insight to the unsettled issue on the T — /3 relation, we have again made 
calculations with the use of methods A and B, by changing M, which is supplementary to 
Ref. y. 

The paper is organized as follows. After briefly reviewing the NES for the quantum 
system, we have derived, in Sec. 2, expressions for the specific heat and susceptibility of 
canonical ensembles of a cluster including dimers both in the BGS and NES, their expressions 
for grand-canonical ensembles having been given in Ref. sl. In Sec. 3 numerical calculations 
of thermodynamical quantities are reported for various values of q and M. The final Sec. 4 
is devoted to discussion and conclusion. In the Appendix, the NES has been applied to a 
cluster of Heisenberg dimers. 

II. NONEXTENSIVE STATISTICS OF HUBBARD DIMERS 
A. Entropy and energy 

We have adopted nanoclusters including independent M Hubbard dimers with negligible 
interdimer interactions, each dimer being described by the two-site Hubbard model. The 
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total Hamiltonian is given by 



M 



H 



(6) 



where denotes the Hamiltonian for the £th dimer, the 
annihilation operator of an electron with spin a on a site j (g £), t the hopping integral, U 
the intraatomic interaction, fis the Bohr magneton and B an applied magnetic field. Six 
eigenvalues of Hi^^ are given by 



a j=l j = l 



[1,2 ei) 



(7) 



ei = 0, 2fiBB, -2i2bB, U, y + A, ^ " A, 



for z = 1 — 6 



(8) 



where A = ^yU^/A + At^ 20| 26|. The number of eigenstates of the total Hamiltonian H is 
First we employ the BGS, in which the canonical partition function for H is given by 



2Q||26| 



Zbg = Tr exp(-/3if). 



^ • ■ ^ exp[-/3(ei, + ■ ■ +eij], 



lyW ]M 
i^BGi ' 



Z'j^l = 1 + 2 cosh{2(3fiBB) + e-""" + 2 e-^^/'cosh(/3A) 



(9) 
(10) 

(11) 
(12) 



where (3 = l/ksT, Tr denotes the trace and Z^^q the partition function for a single dimer. By 
using the standard method in the BGS, we can obtain various thermodynamical quantities 



of the system 



20 



271 ]. Because of the product expression given by Eq. (11), the energy 



and entropy are proportional to M: Ebg = ^-^_bg ^^'^ ^bg = ^S^bg where i?^^ and S^^q 
are for a single dimer. This is not the case in the NES as will be discussed below. 
The entropy Sq in the Tsallis NES is defined by 

m 



Sn 



Tr {pl) - 1 



(13) 
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Here pq stands for the generalized canonical density matrix, whose explicit form will be 
determined shortly [Eq. (16)]. We impose the two constraints given by 



Tr {pm) 



< H >q= Eg, 



(14) 
(15) 



Tr (p^) 

where the normalized formalism is adopted [2|. The variational condition for the entropy 
with the two constraints given by Eqs. (14) and (15) yields 



Pq = ^expg 



with 



Xg = Ti ( exp 



(3 



(3 



{H-Eg 



{H-Eg 



Tr {pl) = Xl-\ 



(16) 

(17) 
(18) 



where exp [x] expresses the g-exponential function defined by 



expjx] = [1 + {l-q)xY- 



0, 



for (1 — q)x > 
otherwises 



and (3 is a. Lagrange multiplier: 



^ dE„ 



(19) 



(20) 



The trace in Eq. (17) or (18) is performed over the 6 eigenvalues, for example, as 



n=l «M=l 



/3 w 



E 



exp„ 



/3 



i^^-Eg) 



(21) 



where the following conventions are adopted: 



i = {k, ■ ■ Hi), 



E 



E--E. 

ii=l «M=1 

H Hi„ 



It is noted that in the limit of g = 1, Eq. (17) reduces to 

X, = Zbg exp[/3Ei] = [ZS exp {^E<^lT. 



(22) 
(23) 
(24) 

(25) 
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For g 7^ 1, however, Xg cannot be expressed as a product form because of the property of 
the g-exponential function: 

expg{x + y) ^ expg{x) expg{y) for g 7^ 1 (26) 

It is necessary to point out that Eg in Eq. (15) includes Xg which is expressed by Eg in 
Eq. (17). Then Eg and Xg have to be determined self-consistently by Eqs. (15)-(19) with 
Eq. (4) or (5) for a given temperature T. The calculation of thermodynamical quantities in 
the NES generally becomes more difficult than that in BGS. 



B. Specific lieat 



The specific heat in the NES is given by 



_ .dl\ (dEg 

" ^ dTj[d/3 



(27) 



Because Eg and Xg are determined by Eqs. (15)-(19), we get simultaneous equations for 
dEq/df5 and dXg/dl3, given by 



with 



ai2 

^21 



X„ 



dEg 
dp 
dXg 



ail 



021 



dEg 
~dj 
dE, 
~d^ 



+ ai2 



dXg 
~dp 



an = g/3Xf2^wf-^e, 



-Xg ^Eg - 



px^ 



9' 



Wi = < i \ exp 



-{-]iH-Eg 



/3 



1 - (1 - g) ^ (e, - i^;, 



l-q 



Wi 



(28) 
(29) 

(30) 

(31) 

(32) 
(33) 
(34) 



(35) 
(36) 
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The specific heat is then given by 



dT / V 1 — ail — fli202i 



with 

-P^, (method A) 

(method B) 



dT 



Xl-'^-l3{l-q)X^\dXjdl3), 
In the hmit of g — )■ 1, Eqs. (27)-(39) yield the specific heat in the BGS, given by [5] 

Cbg = ^ = ksP\< e\ >i - < >\\ 
where < ■ >i is defined by Eq. (15). 



C. Magnetization 



We discuss the field-dependent magnetization, which is given in the BGS, by 

"^BG = ^ = (/^i)l, 

_ sinh(2/3S) 

Zbg 

where fii = —dei/dB, Zbg and < ■ >i are given by Eqs. (9) and (15), respectively. 
In the NES, the magnetization is given by 



dEg ^, ^dSn 

dE„ „_,.._.dX„ 



dB ''OB 

By using Eqs. (15)-(19), we get the simultaneous equations for dEg/dB and dXg/dB 
by 

dE, fdEA^ (^^q\^. 



OB \dB J V dB 

dX. fdEA fdX 

a2l ^TTT +022 ^TTT + 



dB \dB I ' \ dB ' ' 



with 

dl = -Xg^^w1ni + PqX^^^'^^wf'^'^eifXi 

i i 

d2 = PX^g-'j2w!fi,, 
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where Ojj = 1,2) are given by Eqs. (30)-(33). We obtain rUg given by 

-ci2 + r '^-^1 - cii) 



mg = — d2, (49) 

V 1 - Cii - C12C21 / 

= J2 ^1 =< >i ■ (50) 

i 

In the hmit of g — )■ 1, Eqs. (47) and (48) reduce to 

di = - {fii)^ + (3 {eiHi)^ , (51) 
d2 = /3X,{fXi),, (52) 

where < ■ >i is given by Eq. (15) with q = 1. By using Eq. (50), we get mi = niBG, which 
is given by Eq. (42). 

D. Susceptibility 

The high-field susceptibihty in the NES is given by 

UB) = ^. (53) 
The zero- field susceptibility Xq{B = 0) is given by js, 6| 

X, = X,{B = 0) = -Ef + /3-^X-«Xf ), (54) 
where ^ = d^EjdB^ \b=o and ^ = d'^XjdB^ \b=o- With the use of Eqs. (15)-(19), 

(2) (2) 

we get simultaneous equations for Eg and Xq given by 

= ani^f + ai2Xf + A, (55) 

Xf = a^ii^f + a22Xf + /2, (56) 

with 

h = -2/3gX^^-2 5^t.rV?, (57) 

i 

h = /3^gXf-)5^^rV?, (58) 

i 

where Ojj = 1,2) are given by Eqs. (30)-(33). From Eqs. (54)-(58), we get 

X, = — = PqXr' E V.' |b=o . (59) 
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In the limit of q — 1, Eq. (59) yields the susceptibility in BGS: 



Xbg = /3 < 



S=0>1, 



(60) 
(61) 




) 



3 + e-^^ + 2e-/5^/2 cosh(/3A) 




III. CALCULATED RESULTS 
A. The q dependence 

We have performed numerical calculations by changing the index q or the size of a cluster 
M in the NES. Simultaneous equations for Eg and Xg given by Eqs. (15)-(19) have been 
solved by using the Newton-Raphson method with initial values of Ei and Xi obtained from 
BGS {q — 1). The magnetic field B in Eq. (8) is set zero in calculating the entropy, energy 
and specific heat. The calculated quantities are given per dimer. 



First we treat the entropic index g as a free parameter for a single dimer. Figures 1(a)- 
1(f) show the temperature dependence of the energy Eg calculated for = 0. Bold solid 
curves in Fig. 1(a), 1(b) and 1(c) show Ei in the BGS calculated for = 0, 5 and 10, 
respectively The ground-state energy at T = is Ei/t = -2.0, -0.70156 and -0.38516 for 
U/t — 0, 5 and 10, respectively. With increasing q value above unity, the gradient of Eq 
is much decreased in the method A, as shown in Figs. l(a)-l(c). This trend is, however, 
much reduced in the method B, as Figs. l(d)-l(f) show. This behavior is more clearly seen 
in the temperature dependence of the specific heat Cg, as will be discussed shortly [Figs. 



Temperature dependences of the entropy for B — are plotted in Figs. 2(a)-2(f). Figures 
2(a), 2(b) and 2(c) express Sq for U/t — 0, 5 and 10, respectively, calculated by the method 
A, and Figs. 2(d)-2(f) those calculated by the method B. Bold curves denote the results for 



3(a)-3(f)]. 



Fig.2 
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the BGS, where the entropy is quickly increased at low temperature when the interaction 
is increased. When the q value is more increased above unity, 5*^ is more rapidly increased 
at very low temperatures and its saturation value at higher temperatures becomes smaller. 
This behavior is commonly realized in the results calculated by the methods A and B. A 
difference between the two results is ostensibly small because Sq shows a saturation at low 
temperatures. 



Fig.?> 



Figures 3(a)-3(f) show the specific heat calculated for 5 = 0. Ci in BGS for U /t = 
shown by the bold solid curve in Fig. 3(a), has a peak at ksT/t ~ 0.65. Figure 3(c) shows 



that for U /t = 10, this peak splits into two. A lower peak arises from low-lying co 



lective 



spin- wave-like excitations while higher one from single-particle excitations [20|, l27|]. For 
intermediate U/t = 5 these two peaks overlap [Fig. 3(b)]. The temperature dependences 
of the specific heat Cq calculated with the use of the method A for U/t = 0, 5 and 10 are 
plotted in Figs. 3(a), 3(b) and 3(c), respectively. We note that when q is larger than unity, 
peaks become broader. Figures 3(d), 3(e) and 3(f) show the temperatures dependence of the 
specific heat Cq calculated by the method B for U/t = 0, 5 and 10, respectively. Although 
general property of the q dependence of the specific heat of the method B is similar to that 
of the method A, the effect of the nonextensivity in the method B becomes smaller than 
that in the method A. 



FigA 



The BGS susceptibility for U/t = has a peak at ksT/t ~ 0.65 as Fig. 4(a) shows. 
With increasing U/t, the magnitude of xbg is enhanced by the interaction, and its peak 



position becomes lower 20|, l27|, as Figs. 4(b) and 4(c) show: the horizontal scale of Fig. 



4(c) is enlarged compared with Figs. 4(a) and 4(b). The temperature dependences of the 
susceptibility Xg calculated by the method A for f//t = 0, 5 and 10 are plotted in Figs. 4(a), 
4(b) and 4(c), respectively. We note that as increasing q above unity, the peak in Xq becomes 
broader. Figures 4(d), 4(e) and 4(f) show the temperature dependence of the susceptibility 
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Xq calculated by the method B for C//i = 0, 5 and 10, respectively. Again the effect of the 
nonextensivity in the method B becomes smaller than that in the method A. 



B. The M dependence 



Fig.5 



In order to study how thermodynamical quantities of a cluster with Hubbard dimers 
depend on its size M, we have made NES calculations, assuming the q value for a given M 
value with the M — q relation given by Eq. (3). Results for M = oo correspond to those 
of the BGS (g = 1). Figures 5(a)-5(d) show the results for non-interacting case olU/t = 0. 
The specific heat and susceptibility shown in Figs. 5(a) and 5(b), have been calculated by 
the method A with q = 2.0, 1.5, and 1.333 for M = 1, 2 and 3, respectively. Figures 5(c) and 
5(d) express Cq and x^, respectively, calculated by the method B. We note that physical 
quantities in a small cluster with M ~ 1 — 3 are rather different from those of bulk-like 
systems with M = oo, although properties of clusters gradually approach those of bulk with 
increasing M. 



Fig.6 



Similar results for finite interaction of U/t = 5 are shown in Figs. 6(a)-6(d). The specific 
heat and susceptibility plotted in Figs. 6(a) and 6(b), respectively, have been calculated by 
the method A for M = 1, 2, 3 and oo. Figures 6(c) and 6(d) show Cq and Xq, respectively, 
calculated by the method B ior U/t — 5. The results for small M are very different from 
those for M = oo. We note that the M dependences of Cq and Xq of Hubbard dimers shown 
in Figs. 6(a)-6(d) are similar to those of spin dimers described by the Heisenberg model 
[Figs. 13 (a)- 13(d)], details being discussed in the Appendix A. This is not surprising because 
the Hubbard model with the strong coupling and the half-filled electron occupancy, reduces 
to the Heisenberg model with the antiferromagnetic exchange interaction. 



Fig.7 
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We have calculated the M dependence of the maximum values of C* and x* and corre- 
sponding temperatures of and T*. Figure 7(a) shows and T*, and Fig. 7(b) depicts 
C* and X*, all of which are plotted against 1/M: solid and dashed lines denote results cal- 
culated by the methods A and B, respectively. It is shown in Fig. 7(a) that with increasing 
1/M, T* calculated by the method A is much increased than that calculated by the method 
B. We note also that with increasing 1/M, Tq of the method B is increased while that of 
the method A is decreased. Figure 7(b) shows that C* in the method A is smaller than that 
in the method B, whereas x*q the method A is the same as that in the method B. 

IV. DISCUSSIONS AND CONCLUSIONS 

Fig_^ 

Although we have discussed the temperature dependence of physical quantities in the 
preceding section, it is worthwhile to study their magnetic-field dependence. Figures 8(a), 
8(b) and 8(c) show the magnetization rUq as a function of the magnetic field B ior U/t = 0, 
1 and 10, respectively, at ksT/t = 1.0 calculated by the method A: results calculated by the 
method B is not so different from them [Figs. 10(a)]. When q is increased above unity, the 
magnetization at lower fields {h/t < 1) is decreased whereas at higher fields {h/t > 1) it is 
much increased. This is consistent with the calculation of the susceptibility shown in Fig. 
3(a)-3(c), where Xq at ksT/t — 1.0 is smaller for larger q: the susceptibility stands for the 
initial gradient of rUq at S = 0. 





Fig.lQ 

The field dependence of physical quantities are discussed in more details. Figure 9 shows 
the B dependence of the six eigenvalues of for C//i = 5 [Eq. (8)]. We note the crossing of 
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the eigenvalues of es and ee at the critical filed: 



4' 



(62) 



leading to ^BBc/t = 0.351 for U/t = 5. At B = the magnetization nig is rapidly increased 
as shown in Figs. 10(a) and 10(b) for ksT/t = 1.0 and 0.1, respectively: the transition at 
lower temperature is more evident than at higher temperature. This level crossing also 
yields a peak in Xq [Figs. 10(c) and 10(d)] and a dip in Cg [Figs. 10(e) and 10(f). It is 
interesting that the peak of Xq for q = 1.5 in the NES is more significant than that in the BGS 
whereas that of Cg of the former is broader than that of the latter. When the temperature 
becomes higher, these peak structures become less evident. Similar phenomenon in the field- 
dependent specific heat and susceptibility have been pointed out in the Heisengerg model 
within the BGS 3- 

Figure 10(a) and 10(b) remind us the quantum tunneling of magnetization observed in 
magnetic molecular clusters such as Mn4, Mnl2 and Fe8 29|. It originates from the level 
crossing of magnetic molecules which are parallel and anti-parallel to the easy axis when a 
magnetic field is applied. 

The N — q relation given byg = l-|-2/A^ [Eq. (2)] has been derived from the average of 
;he BGS partition function of exp(— /3e) with e = J2i ^« o'^^r fluctuating (3 fields, as given by 
2ll-|23| 



N 

d(3 exp ( — /3 e,; 

1=1 



= exp, 



N 



i=l 



with 



r(/3) 

2 

N 



r(f) Wo) ^ 

</3>, 

</3^> - <(3>^ 



exp 



Wo 



(63) 

(64) 
(65) 
(66) 



< /3 >2 

Here < Q > stands for the expectation value of Q averaged over the F (or x"^) distribution 
function f^{P), the average of the fluctuating /3 and 2/N its variance. The F distribution 
is emerging from the sum of squares of Gaussian random variables. Alternatively, by using 
the large-deviation approximation, Touchette 30|| has obtained the distribution function 
/'^(/3), in place of /^(/3), given by 



/3o fN(3o 
F(f) 



/3- 



exp 



m 

' 2/3 



(67) 
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Fig.ll 



Figure 11 shows the f^- and /^-distribution functions for various values. For N oo, 
both reduce to the delta-function densities, and for a large = 100, both distribution 
functions lead to similar results. For a small (< 10), however, there is a clear difference 
between the two distribution functions. We note that a change of variable /3 — )■ in /'^ 
yields the distribution function similar to It should be noted that /"^ cannot lead to the 
g-exponential function which plays a crucial role in the NES. For a large e, leads to the 
stretched exponential form of w{t) ~ e'^^ while yields the power form of w{e) ~ e^^^. 
This issue of vs. is related to the super statistics, which is currently studied with much 
interest |c 

In deriving Eq. (63), we have implicitly assumed that the distance between sparsely 
distributed clusters is larger than ^, the coherence length of the fluctuating (5 field over 
which the field /3 uniformly fluctuates, and that the linear size of the clusters is smaller than 
^. If the population of constituting dimers in a given cluster is sparse such that the distance 
between dimers is larger than ^ and a local fluctuating /3j field around a dimer i is almost 
independent from the local /3j field around another dimer j, Eq. (64) is replaced by 

w({e,}) = nil, / rfA exp(-AQ) = niI,exp^,[-/3oe,], (68) 

with g' = 3. Actual distribution of dimers is expected to lie between the two extreme cases 
given by Eqs. (63) and (68) [22|. 

Although results calculated by the two methods A and B are qualitatively similar, there 
are some quantitative difference, as previously obtained in Refs. sl and js]. When we 
calculate the Curie constant F^ of the susceptibility defined by Xq{T) = {fi%/kB)[Tq{T)/T], 
the ratio between F^^^ and Tq^^ calculated by the two methods, is given by 

Jb) = (69) 
= 4^''-^), for T = (70) 

= 6^"^^). for T = oo (71) 

In general, Tg^^ depends on t, U and T. In the limit of t = 0, for example, it is given by 32 1 

Ff ) = 2 g, for T = (72) 

4 

= - g. for T = oo (73) 
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The result of the method A given by Eqs. (69)-(73) leads to anomalously large Curie constant 
compared to that of the method B. Equation (69) is consistent with^the result for spin dimers 
described by the Heisenberg model (for details, see Appendix A) [5]. A comparison between 
Eqs. (35) and (63) yield the average temperature < T > given by 

^ /3o = ^ (74) 



kB<T> (3 

which supports the method B. These results suggest that the method B is more appropriate 



331, 



3^ (for 



than the method A. This is consistent with recent theoretical analyses made in 
related discussions, see Ref. |35|]. 

In summary, within the framework of the NES, thermodynamical properties have been 
discussed of a cluster including M dimers, each of which is described by the two-site Hubbard 
model. We have demonstrated that the thermodynamical properties of small-scale systems 
are rather different from those of bulk systems. Owing to recent progress in atomic engineer- 
ing, it is possible to synthesize molecules containing relatively small numbers of magnetic 
atoms with the use of various methods (for reviews, see refs. 8-10). Small-size magnetic 
systems ranging from grains (micros), nanosystems, molecular magnets and atomic clusters, 
display a variety of intrigue physical properties. It is interesting to compare our theoretical 
prediction with experimental results for samples containing small number of transition-metal 
dimers of M = 1, 2 and 3. Unfortunately experiments on samples with such a very small 
number of dimers, have not been reported. Theoretical and experimental studies on nan- 
oclusters with changing M could clarify a link between the behavior of the low- dimensional 
infinite systems and nanoscale finite-size systems. The unsettled issues on T — /3 and the 
N — q relations in the current NES are expected to be resolved by future theoretical and 
experimental studies on nanosystems, which are expected to be one of ideal systems for a 
study on the NES. Our discussion in this study has been confined to the static property 
of nanoclusters. It would be interesting to investigate dynamics of dimers which has been 
discussed within the framework of the BGS. 
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Appendix: NES for Heisenberg dimers 



We consider a cluster containing M Heisenberg dimers, spin dimers described by the 
Heisenberg model (s = 1/2), given by 

M 

H = Y,Hf\ (Al) 



=1 



Hf = - Jsi ■ S2 - giJiBB{s^, + S2.), (1, 2 e £) (A2) 

where J stands for the exchange interaction, g (=2) the g-factor, the Bohr magneton, 
and B an applied magnetic field. Four eigenvalues for are given by 

ei = — g^BBrrti, with mi = 1, 0, —1 for z = 1, 2, 3 

3J 

= — — giisBmi. with rrii = for i = 4 (A3) 

The number of eigenvalues of H becomes 4*^. 

In the BGS the canonical partition function is given by 

Zbg = [Z^B%r, (A4) 
= exp [1 + 2cosh {2/3fiBB)] + exp , (A5) 

with which thermodynamical quantities are easily calculated. The susceptibility is, for 



example, given by 



12 



13| 



XBG = Mxfc, (A6) 

""^'^ ^ {k^) {s + expi-J/kBT)) ■ ^^^^ 

The calculation of thermodynamical quantities in the NES for the Heisenberg model goes 
parallel to that discussed in Sec. 2 if we employ eigenvalues given by Eq. (A3). For example, 
by using Eq. (59), we get the zero-field susceptibility for Heisenberg dimers, given by 

X, = 9'f^l(f) ^j:^r'ml (A8) 

\Cg / u\q . 

where a sum Y2i is performed over 4*^ eigenvalues [see Eq. (23)]. In the case of M = 1 (a 
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single dimer), we get 



2 2 I'^lP 



exp. 



/3 



J 



2g-l 



3exp„ 



+ exp 



(A9) 
(AlO) 



3J\ 



exp„ 



}• 



(All) 



In the limit of g = 1, Eq. (A9) reduces to given by Eq. (A7). 

The Curie constant Eg defined by Xq = {I-'''bI^b){^ q/T) for J = is given by jsl 



2Mg, 



(method A) 
(method B) 



which are consistent with results obtained for Hubbard dimes [S]. Equations (A12) 
an anomalously large Curie constant, which was referred to as dark magnetism in Ref 



(A12) 
(A13) 

leads to 



Fig.l2 



Figure 12(a) and 12(b) show the temperature dependences of the specific heat and suscep- 
tibility, respectively, of a single Heisenberg dimer (M = 1) for several q values for J < (an- 
tiferromagnetic coupling) calculated by the method A (solid curves) and B (dashed curves). 
Figures 12(c) and 12(d) show the specific heat and susceptibility, respectively, of two spin 
dimers (M = 2) calculated with the use of 4^ eigenvalues. 



Fig.l?> 



Figures 13(a) and 13(b) show Cq and Xq when the size M of a cluster of Heisenberg 
dimers is changed, with q = 2.0, 1.5, 1.333 and 1.25 for M = 1, 2, 3 and 4, respectively, 
calculated by the method A: results for g = 1 of the BGS [corresponding to M = oo in Eq. 
(3)] are included for a comparison. Figs. 13(c) and 13(d) show similar results of Cq and Xq 
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calculated by the method B. The M dependence of Cq and Xq fo^^ Heisenberg dimers shown 
in Fig. 13(a)-13(d) is quite similar to those shown in Figs. 6(a)-6(d) for Hubbard dimers. 
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FIG. 1: The temperature dependences of the energy Eq for B = with (a) U/t = 0, (b) 5 and 
(c) 10 calculated by the method A, and those with (d) U/t = 0, (e) 5 and (f) 10 calculated by the 
method B: q = 1.0 (bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) 
and 2.0 (solid curves). 

FIG. 2: The temperature dependences of the entropy Sq for B = with (a) U/t = 0, (b) 5 and 
(c) 10 calculated by the method A, and those with (d) U/t = 0, (e) 5 and (f) 10 calculated by the 
method B: q = 1.0 (bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) 
and 2.0 (solid curves). 

FIG. 4: The temperature dependences of the susceptibility Xq for S = for (a) U/t = 0, (b) 5 and 

(c) 10 calculated by the method A, and those for (d) U/t = 0, (c) 5 and (f) 10 calculated by the 
method B: q = 1.0 (bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) 
and 2.0 (solid curves). 

FIG. 5: The temperature dependences of (a) specific heat Cq and (b) susceptibility Xq (per dimer) 
of Hubbard dimers for U/t = calculated by the method A, and those of (c) specific heat Cq and 

(d) susceptibility Xq calculated by the method B, with M = 1 (bold solid curves), M = 2 (chain 
curves), M = 3 (dashed curves) and M = oo (solid curves). 

FIG. 6: The temperature dependences of (a) specific heat Cq and (b) susceptibility Xq (per dimer) 
of Hubbard dimers for U/t = 5 calculated by the method A, and those of (c) specific heat Cq and 
(d) susceptibility Xq calculated by the method B, with M = 1 (bold solid curves), M = 2 (chain 
curves), M = 3 (dashed curves) and M = oo (solid curves). 

FIG. 3: The temperature dependences of the specific heat Cq for S = with (a) U/t = 0, (b) 5 
and (c) 10 calculated by the method A, and those with (d) U/t = 0, (e) 5 and (f) 10 calculated 
by the method B: q = 1.0 (bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain 
curves) and 2.0 (solid curves). 
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FIG. 7: (a) 1 jM dependence of the temperatures of (circles) and T* (squares) where Cg and 
Xq have the maximum values, respectively, (b) 1/M dependence of the maximum values of C* 
(circles) and %* (squares). Solid and dashed lines denote the results calculated by the methods A 
and B, respectively: T* calculated by the method A shown in (a) is divided by a factor of five. 



FIG. 8: The magnetization as a function of the magnetic field B for (a) U/t = (b) 1 and 
(c) 10 at ksT/t = 1.0: q = 0.8 (double-chain curves), 0.9 (chain curves), 1.0 (bold solid curves), 
1.1 (dotted curves), 1.2 (dashed curves), and 1.5 (solid curves) calculated by the method A. 



FIG. 9: The magnetic-filed dependence of the eigenvalues Cj (i = 1 — 6), for U/t=5. 



FIG. 10: The magnetic-filed dependence of (a) the magnetization rriq for ksT/t = 1.0 and (b) 
ksT/t = 0.1, (c) the susceptibility for ksT/t = 1.0 and (d) ksT/t = 0.1, (e) the specific heat Xq 
for ksT/t = 1.0 and (f) ksT/t = 0.1, of a single Hubbard dimer (M = 1) for U/t=5, calculated by 
the method A (solid curves) method B in the NES (dashed curves), and in the BGS (chain curves). 



FIG. 11: The distributions of f^{P) (solid curves) and /'^(/3) (dashed curves) as a function of P 
(see text). 
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